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I ■ Abstract 

Oh! 

D . A simple model for the instanton ensemble at finite temperature T is proposed, 



including "random" and strongly correlated "molecular" component. T-dependence 



S^ . of fermionic zero modes naturally leads to chiral symmetry restoration, without 

C^ \ instanton suppression. Moreover, at T = (1 — 2)Tc the non-perturbative effects due 

to "molecules" are so strong, that they even dominate the global thermodynamics. 
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Instantons are the major component of non-perturbative fields in the QCD vac- 
uum, and significant amount of work |||, ^, |3| has been done in order to build a 
quantitative theory describing them. Two major steps forward were done during 
the last year. First, the simplest ensemble (the so called Random Instanton Liquid 
Model), has reproduced many correlation functions Q known from phenomenology 
[^ and lattice simulations p. Second, by "cooling" of the (quenched) lattice config- 
urations it was found in [^ that the typical instanton density is about n ~ lAfm~'^ 
and the typical size is about p ~ .35/m, very close to the "instanton liquid" pa- 
rameters suggested by one of us a decade ago 0. The correlation functions and 
hadronic wave functions, are shown to be practically unaffected by "cooling": so, 
by removing perturbative gluons and confinement, one does not loose mesons and 
baryons! 

The instanton-induced effects at finite temperature T and much less studied. 
Our particular focus is on the chiral restoration phase transition at Tc- The main 
idea is that it happens due to growing correlations between instantons and anti- 
instantons, described by the disappearance of "single" instantons and growth of // 
"molecules" . 

Our first study along this line was performed few years ago ||8[. However, the 
present paper is significantly different, because it is based on completely new mech- 
anism. Instead of Debye-type screening P|, implemented by thermal suppression 
factor f{T) = exp{— p^T"^ const) [|lO| Q , the phase transition occurs due to the 
T-dependence of the quark-induced II interaction. 

The ensemble of interacting instantons can be studied with the partition func- 
^This suppression is expected to work at high T » Tc only, and available lattice data on screen- 
ing mass suggest that it may probably happens at only at temperatures T > 300 — 400Mel/. Direct 
lattice measurements of the instanton density |ll| have indeed seen no significant T dependence 
till such high temperatures. 



tion 

N++N. 

^=/ E Wl^n n dn.d{p^)p^^exp{-S,nt)Il^'det{iD + imf) (1) 

where dQi is the measure in space of cohective coordinates, (12 per instanton in 
SU{3)), d{p) the instanton amphtude, and Smt the gluonic interaction. The last 
factor, appearing after integration over fermions, is the one we are focussing on. 
Assuming A^+ = N^ = N, one writes it in terms of a iV x A^ 'hopping' matrix, 
with TjA = J d'^x{ipAo{x ~ ^a) \iDx {tpioix — zj)) from some instanton / to some anti- 
instanton A {za, zi are the centers of j4, /, ^o the zero-modes), as det{TT^ + m'j,). 
The statistical system described by this partition function is quite complicated, and 
direct simulations were done so far only for T = 0. 

The instanton solutions and their zero modes are known analytically for non- 
zero T, and detailed studies of 'hopping' matrix elements were done in |12|. They 



have the structure Tja = 'U4/1 + (ur/r)/2, where the 2x2 matrix u^t^ describes 
the relative orientation of / and A. Rather complicated formulae for /i,/2 were 
derived in [12|. 



We consider the instanton ensemble as a superposition of the uncorrelated (or 
"random"') component and a highly correlated (or "molecular") component, with 
the (4-d) densities na{T) , nra{T) . In the former component instantons are assumed 
to have random relative orientation u^, and, as in |p, it is treated in a mean field 
approximation. The fermionic determinant is obtained from /(T) -^ (TjaT^j), i.e. 
summed over anti-instantons A with a density na/2, with random positions and 
orientations. 

For the "molecular" component we assume the opposite, namely the most fa- 
vorable relative color orientation u^ ~ {zj — z^), maximizing the hopping matrix 
element \Tja\ (as well as exp{—Sint))- Therefore, we have to calculate another func- 
tion I(Nf,T) ~ {{TiaT\i) 1) where now averaging means integration only over the 
relative coordinates inside a II pair. 



Temperature dependence of these average matrix elements was obtained by nu- 



merical integration of the formulae from |12], and the results are shown in Fig. 1. 



One has to exclude too close II pairsg , and we show results with the "core radii" 
Re = 1/9 (dashed) and Re = 2p (solid). Although the overlap integrals significantly 
depend on Re, the resulting uncertainty of the thermodynamical quantities (see be- 
low) is in fact not so dramatic. Note remarkably different T-dependence of these 
two quantities. While the I{T) (indicated as the Nj = 1 curves in Fig.l) decreases 
with T, I{Nf,T) grows, the stronger the larger Nf. As a result, the molecular com- 
ponent builds up with T, while the random one decreases and eventually disappears, 
restoring the chiral symmetry. 

In order to simplify gluonic interaction, and still describe the system self- 
consistently, we adopt a simple "average repulsion" (Sint) = i^P^P^ with one di- 
mensionless parameter k. The same parametrization is used for the description 
of interaction between all instantons, belonging to the random component or to 
molecules. 

Let us now evaluate the statistical sum of the system in terms of the densities 
na,nm, starting with the differential activities for molecular 

dZm = C^dpidp2<fR'^^^^^{piP2f-''exp[-K{pl + pl){plna + 2/5^71^)] (T/^Tl,)^/ 
i ^5(7(3) 

and random components 

dza = 2Cdpp^-''exp[-Kp^{plna + 2plnm)] < TT+ >^f (3) 

where b = -^-Nc — nNj is the coefficient of the Gell-Mann-Low function and 

< TT^ >= p^'^[i/(T) / dna{p)p\^''^ . Besides the densities, an important ingredient 

of the interaction are the root mean square radii pa,Pm, which can be found from 



^The general reason for such 'repulsive core' is generally related to the fact, that too close pairs 
do not in fact correspond to strong fields and are not objects of the semiclassical theory. 



eqs. (3) and (4) to be related to each other through 

4 = ^; a = 6/2-l; /? = 6/2 + 3iV^/4 - 2 (4) 

Pa P 

Another relation between them connects the interaction parameter k to the dilute- 
ness of the ensemble 

1 _ 2p^na 4p^ra^ 

^"^r + ^^ • ^^> 

So, one can eliminate the mean square radii and get the activities 



[n„ + (2a//3)n„]°' (4k/3)" 



" [n„ + (2a//3)n„]^/2+^//8' (2ac)/^ ^ 2 ^ ^k^ 

As usual, the grand potential 



(7) 



n = -p= -{logZ)/V, = ^logC-^) + ^logC-^) (8) 

(where V4 is the 4-dimensional volume) should then be minimized with respect to 
the particle numbers Na = V^na and Nm = V/j^Um. The resulting grand potential 
provides the instanton contribution to the pressure p and to the energy density 
e = -p + T^. 

The general case leads to rather cumbersome resulting equations, but for T > Tc 
one has Na = and the following simple result p = 2^™. 0- 

The statistical sum under consideration describes two phases, with and without 
chiral symmetry. However, the realistic description should include the contributions 
to thermodynamics unrelated with instantons as well: it is the total pressure, which 



^ It is interesting, that connecting Um with the gluon condensate ^^(GJJ^GJi^) — Wrim one 
reproduces the famous "trace anomaly" expression. However, that the contribution of molecules 
to the energy density is not just the same expression with the opposite sign, because the fermionic 
determinant has an explicit T-dependence. 



should be continuous through the transition. We use the simplest possible model 
here, including the non-interacting massless pions in the broken phase (which for 
Nf = 2,3 are actually irrelevant) and the ideal quark-gluon plasma in the symmetric 
phase. 

A typical set of results is shown in Fig. 2 for Nf = 2 and the cores ^ = I and 2, 
to see uncertainties involved. The constants A and B entering the activities could be 
determined from first principles, provided we know the accurate value of Aqqd and 
instanton interactions. We select B in order to get the instanton density at T = 
equal to that found on the lattice, lAfm~'^, and A in order to get Tc — 150MeV. 

The upper panel of Fig. 2 tells us, that although the "random component" 
(solid line) dominates the broken phase, the number of molecules jumps up at 
the transition, and the total instanton density above Tc, 2njn, turns out to be 
comparable to that at T = 0. We have not plotted the behaviour of the quark 
condensate, which scales as < tpip >'^ UaiTY''^ . Inside uncertainties of the model 
it is essentially flat, till nearly T^, which qualitatively agree with lattice data. 

The importance of the molecular component is better demonstrated by the T- 
dependence of the pressure (the middle panel of Fig. 2): in fact, for T = (1...2)Tc 
the contribution of "molecules" (dash-dotted line) is crucial, without it one would 
not be able to sustain the pressure because that from quarks and gluons (the dotted 
line) is not large enough ! 

The lower panel shows the large jump in energy density at the transition. Al- 
though it is mostly due to the "liberation" of quarks and gluons, a finite part of it 
is also generated by the "molecules" (dash-dotted line). Such behaviour was in fact 
observed in lattice calculations with dynamical fermions: unlike those for the pure 
glue case, they show e(T) which overshoots the Stefan-Boltzmann value above Tc- 

Finally, let us discuss what happens at larger Nf. First of all, from (6) it is 
^ Note that both pressure and energy density are counted not from physical, but from pertur- 
bative vacuum. 



clear that a positive radius of molecules can only be obtained for a > 0, 6 > 2, or 
Nf < 13.5 [^]. Furthermore, for fixed interaction (namely for core parameter Re 
and K being independent on Nf) the broken phase exists in a shrinking region of 
temperatures, and for Nf > N'^pp'^'^ the model predicts the unbroken chiral symme- 
try even in the ground state of the theory. The vacuum structure for large Nf was 
studied on the lattice, and (although the question is by no means settled) indica- 
tions were reported |l^ that N^^^^"^ = 7. This number is similar to what one gets 
in our model. 
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Figure Captions 

Fig. 1. The ratio I{Nf,T)/I{Nf,0) of the average fermionic determinant for 
molecules at temperature T, normalized to its value at T=0. Numbers on the figure 
are Nf, the number of (massless) quark flavors. The curves marked by 1 correspond 
to /(T)//(0), the integral for "random" component. Solid curves are for core radius 
Re = 2/9, and the dashed ones are for Re = p. 

Fig. 2. Chiral restoration phase transitions for 2 massless flavors and two differ- 
ent core parameters Re- The upper panel shows the T dependence of the densities Ha 
(solid) and Um (dotted). In the middle panel the T dependence of the total pressure 
p (solid) is shown, including the contributions of the pion gas/quark-gluon plasma 
(dotted) and of instantons (dash-dotted). The energy density is presented in the 
lower panel (solid) which is modified by the instanton contribution (dash-dotted). 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9401285vl 



This figure "fig2-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9401285vl 



